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SHEAR  EFFECTS  ON  OCEAN  ACOUSTIC  PROPAGATION  DUE  TO  STEP-PERIODIC 
ROUGHNESS  ALONG  THE  OCEAN  BOTTOM  INTERFACE 

Stanley  A.  Chin-Bing  t  and  Joseph  E.  Murphy  § 

t  Numerical  Modeling  Division,  Naval  Oceanographic  and  Atmospheric 
Research  Laboratory,  Stennis  Space  Center,  MS  39529-5004,  U.S.A. 

§  Department  of  Physics,  University  of  New  Orleans,  New  Orleans,  LA  70148,  U.S.A. 


ABSTRACT 

Low-frequency  ocean  acoustic  propagation  can  be  affected  by  backscatter  from  rough  ocean 
bottoms  and  by  compressional-shear  wave  conversions  within  the  ocean  bottom.  Numerical 
simulations  have  been  performed  that  include  the  effects  of  such  backscattering  and  loss 
mechanisms.  A  full  wave  range-dependent  finite-element  ocean  seismo-acoustic  computer  model 
(SAFE)  was  used  together  with  a  simulated  rough  ocean  bottom  structure  that  supported  compres- 
sional  and  shear  wave  speeds  as  well  as  compressional  and  shear  wave  attenuations.  Step-periodic 
changes  in  depth  and  range  along  the  water-bottom  interface  were  used  to  simulate  ocean-bottom 
roughness.  A  numerical  example  was  used  to  illustrate  the  combined  effects  of  ocean-bottom 
shear  and  roughness,  and  to  give  insight  into  these  complicated  processes. 

KEYWORDS 

Shear  waves;  seismo-acoustics;  ocean  acoustics;  wave  propagation;  backscatter. 

INTRODUCTION 

Recently,  the  Finite  element  (FE)  method  has  been  applied  (Murphy  and  Chin-Bing,  1988.  1989, 
1991;  Chin-Bing  and  Murphy,  1988,  1990)  to  develop  several  computer  models  that  calculate  the 
complex  pressure  field  resulting  from  acoustic  and  elastic  propagation  and  scattering  in  a 
inhomogeneous  ocean  with  range-varying  boundaries  (sea  surface  and  ocean  bottom).  Although 
these  computer  models  are  computational  intensive,  they  have  been  verified  to  give  benchmark 
accurate  predictions  of  the  underwater  complex  pressure  field.  In  this  paper  we  first  briefly  review 
the  salient  points  leading  to  the  developments  of  these  FE  acoustic  and  elastic  ocean  propagation 
models.  Next,  further  developments  and  implementations  regarding  our  seismo-acoustic  FE 
model.  SAFE  (Murphy  and  Chin-Bing,  1991),  are  presented.  The  SAFE  model  includes  the  ocean 
bottom  as  an  aneiastic  medium,  having  range  and  depth  dependence  in  compressional  and  shear 
speeds  and  in  compressional  and  shear  attenuations.  Finally,  we  apply  the  SAFE  model  to 
simulate  a  typical  shallow  water  scenario  with  a  rough  ocean-bottom  interface  overlying  a  shear 
supporting  media.  Further  verification  of  the  accuracy  of  the  SAFE  model  is  also  presented. 

In  the  development  of  the  scalar  FE  ocean  acoustic  model,  FOAM  (Murphy  and  Chin-Bing.  1988. 
1989),  the  starting  point  is  the  inhomogeneous  Helmholtz,  equation, 

pV*  (p-1  VP)  +  k2  P  =  -  6(r  -  rs)  .  (1) 

The  solution  to  this  equation  (with  appropriate  boundary  conditions)  gives  the  time-harmonic 
complex  pressure  field  from  a  unit-strength  point-source  of  angular  frequency  a)  at  position  rs  in 
an  inhomogeneous  ocean.  A  three-dimensional  ocean  environment,  axially  symmetric  about  a 
vertical  line  passing  through  the  source,  is  assumed,  In  Eq.  (1),  written  in  cylindrical  coordinates, 
p  »  p(z,  r)  is  the  density,  k  -  k(z,  r)  ■  to/c(z,  r)  is  the  wavenumber,  and  c(z,  r)  is  the  sound  speed. 
Attenuation  is  included  as  an  imaginary  term  in  the  complex  wavenumber. 

FOAM  simulates  a  typical  ocean  acoustic  scenario  as  follows:  The  acoustic  source  is  arbitrarily 
placed  in  the  water  column  or  fluid  ocean  bottom;  the  sound  speeds,  in  the  water  and  in  the  fluid 


780 


Proc .  of  the  8th  ICMCM 


7*1 


ocean  bottom,  may  vary  with  both  depth  and  range;  the  depth  of  the  water-bottom  interface  may 
vary  with  range;  the  acoustic  pressure  vanishes  at  the  air-water  interface  (pressure  release 
surtace);  beyond  some  maximum  depth  the  ocean  bottom  is  assumed  to  be  homogeneous;  beyond 
some  maximum  range  the  sound  speed  profiles  are  assumed  not  to  vary;  then,  along  the  right  and 
bottom  computational  boundaries,  FOAM  imposes  an  approximate  boundary  condition  (radiation 
condition)  to  prevent  energy  reflecting  from  the  computational  boundaries  back  into  the  interior. 

In  using  the  FE  method  to  solving  the  ocean  acoustic  (fluid)  wave  equation,  the  differential 
equation  and  its  unknown  solution  are  replaced  by  a  system  of  algebraic  equations  in  terms  of 
parameters  that  define  an  approximate  solution.  The  domain  is  partitioned  into  non-overlapping 
elements  and  assumes  a  simple  form  for  the  approximate  solution  within  each  element.  These 
local  representations  are  then  joined  together  by  using  the  appropriate  physical  continuity 
conditions  to  provide  a  global  solution.  The  ocean  acoustic  wave  equation  and  the  associated 
boundary  conditions  give  a  system  of  matrix  equations  which  are  sparse,  banded,  and  symmetric. 

FOAM  uses  a  FE  grid  composed  of  triangular  elements  and  computes  the  complex  pressure 
values  at  the  each  node  of  the  grid.  The  FE  grid  is  formed  so  that  the  boundaries  of  the 
appropriate  interior  elements  coincide  with  the  actual  physical  interfaces,  thus  giving  a  more 
physically  correct  computational  model.  Discontinuities  in  density,  sound  speeds,  etc.  that 
normally  occur  in  going  across  the  water-bottom  interface  are  thus  easily  included  in  FOAM. 


THE  SEISMO-ACOUSTIC  FINITE-ELEMENT  MODEL:  SAFE 


The  FE  methods  employed  in  creating  the  underwater  acoustic  (compressional)  model,  FOAM, 
have  been  used  to  create  the  Seismo-Acoustic  Finite  Element  model,  SAFE,  that  includes  the 
coupled  effects  of  shear  waves  and  compressional  waves,  together  with  shear  and  compressional 
attenuations.  The  SAFE  model  departs  from  FOAM  in  three  major  ways:  (1)  Rectangular 
elements  are  used  (Fig.  1)  rather  than  the  triangular  elements  used  in  FOAM;  (2)  Displacements 
are  calculated  by  the  SAFE  model  and  the  complex  pressure  field  obtained  from  derivatives  of 
these  displacements  -  FOAM  solves  for  the  complex  pressures  directly;  and,  (3)  SAFE  has  an 
inhomogeneous  anelastic  ocean  bottom  -  FOAM  assumes  an  inhomogeneous  fluid  ocean  bottom. 
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Fig.  1.  SAFE  finite  element  grid  with  boundary  conditions.  Element 
nodes  are  adjusted  to  lie  along  the  water-bottom  interface. 


2-D  Elastic  Equations  of  Motion 

The  equations  for  a  2-D  elastic  medium  (in  the  xy-plane)  can  be  written  in  matrix  form  as, 

f  s-a^p  d  -  DTCDd  .  (2) 

t  =  NT  C  D  d  ,  (3) 

where:  d  =  [  u.  v  ]T  is  a  column  matrix  containing  node  displacements  u.  and  v;  f  =■  [  fx,  fy  ]T  is  a 
column  matrix  containing  horizontal  and  vertical  body  forces,  fx,  fy  ,  respectively;  C  is  a  square 
matrix  containing  the  elasticity  constants,  and  t  =  [  tx,  ty  JT  is  a  column  matrix  composed  of  the 
boundary  (traction)  forces.  D  is  a  matrix  of  differential  operators  (partial  derivatives  with  respect 
to  the  x  and  y  coordinates)  and  N  is  a  matrix  composed  of  x  and  y  components  of  the  unit  normal 
vector,  n  =  (  n„,  ny ),  to  the  boundary.  They  are  defined  by: 
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r)  /  Sx  O  ()  /  <h 

(4) 

0  <)  /  th  D  /  <)x 

For  plane  strain.  Cn  =  c22  =  E  (i-v)  /  (  (i+v)(i-2v)  ]  ,  C12  =  Ev  /  [  (1+v)(i-2v)  j  ,  and  c33  = 
Ev  /2(i+v)  E  is  the  modulus  of  elasticity,  and  v  is  Poisson's  ratio. 

The  equations  of  motion  for  plane  elasticity  (equilibrium  of  stresses,  strain-displacement 
relations,  and  stress-strain  or  constitutive  relations)  are  represented  by  Eq.  (2)  and  the  boundary 
conditions  (i.e.,  traction  forces  on  the  boundaries)  are  represented  by  Eq.  (3).  In  the  SAFE  model, 
axial  symmetry  is  assumed  and  Eqs.  (2)  -  (4)  are  written  appropriately  in  cylindrical  coordinates. 

The  bilinear  interpolations  functions  for  the  rectangular  element  are  =  (i-^.a)(l-ri/b), 

T2(^.n)  =  (4/a)(i-r|/b),  =  (^/a)(n/b),  and  T4(£.t|)  =  (l-^/a)(n/b),  where  £,  and  n  arc  local 

coordinates,  £  =  x-xi,n  =  y-yi.  and  T  obeys  the  orthogonality  conditions  and  is  normalized,  i.e., 
T,  (^i/nj )  =  Sjj  and  E  Hq  )  =1,  where  the  summation  extends  from  1  to  4.  Using  the  finite 
element  interpolation  for  u  and  v  in  each  individual  element,  (i.e.,  u  =  L  U;  Tj ,  v  =  z.  vt  lP,  .with 
the  summation  extending  from  1  to  4)  yields  matrix  equations  of  the  form. 


M. 


n. 


I)  = 


(K-  M  )  U  =  F  , 


(5) 


with,  K  =  JnBTCBdxdy,  F  =  J^fdxdy  +|r  T^dS,  U  =  [u„ 


and.  B  =  DVF  , 


M  =  {Dip'Vr'V  ,  = 
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%  0  % 
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ro 

,  v2 
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0  ' 
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(6) 
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Different  material  properties  are  allowed  within  each  rectangular  element.  The  material  constants. 
Cjj ,  and  the  density  are  kept  constant  within  an  element,  thus  allowing  all  of  the  integrations 
required  for  the  K-matrices  to  be  done  analytically.  One  could  include  the  ability  to  deform  the 
rectangular  elements  so  that  the  nodes  could  follow  the  range-dependent  ocean  bottom 
bathymetry,  or  any  other  surfaces  across  which  there  are  sharp  discontinuities  in  properties. 
However,  numerical  integrations  would  then  have  to  be  made  for  tlr  contributions  that  these 
quadrilateral  elements  would  make  to  the  FE  system  of  equations. 

The  Boundary  integral  andJRadiaiLngJ.Guadiiry.Cgndiiions 

The  FE  rectangular  elements  are  influenced  by  neighboring  elements  and  boundaries  depending 
upon  whether  the  elements  are  interior  elements  or  elements  with  one  or  two  of  its  sides  on  a 
boundary.  Elements  within  the  interior  of  the  FE  grid  are  of  the  form,  I  T,  t*  dS  and 
/  T,ty  dS  .  These  integrals  are  exactly  cancelled  by  similar  terms  from  a  neighboring  element. 
This  is  due  to  the  condition  that  the  stresses  must  be  continuous.  Elements  with  one  or  possibly 
two  of  their  sides  on  the  actual  boundary  of  the  computer  model  are  allowed  two  different 
boundary  conditions:  (a)  Stress-Free:  t*  =  ty  =  0.  This  is  treated  exactly  and  the  boundary  term  is 
zero,  (b)  Radiation  Boundary  Condition:  Waves  should  propagate  through  the  model's 
terminating  boundaries  without  reflection.  This  boundary  condition  is  treated  approximately. 

In  the  scalar  CW  pressure  wave  model,  FOAM,  the  normal  derivative,  9P/9n,  appears  in 
boundary  integrals  computed  along  the  bottom  and  the  right  side  of  the  computational  grid.  The 
approximate  radiation  condition  First  used  was  dP/9n  -  ikP  »  i(a>  zc)P,  which  gives  a  small  effective 
reflection  from  the  boundary.  In  the  seismo-acoustic  model,  SAFE,  approximate  boundary 
conditions  discussed  by  Kalinowski  (DiNapoli,  and  Deavenport,  1979)  are  incorporated:  on  = 
p  eg  vn  (for  norma!  traction),  and  cq  =  p  cs  vt  (for  tangential  traction),  where  on  ,  ot  are  the 
normal  and  tangential  stress,  respectively;  Cp,  cs  are  the  compressional  and  the  shear  wave  speeds, 
respectively;  and,  vn  ,  vt  are  the  normal  and  the  tangential  components  of  panicle  velocity  on  the 
boundary,  respectively.  These  conditions  for  both  the  scalar  wave  model  and  the  seismo-acoustic 
model  are  exact  for  plane  waves  propagating  normal  to  the  boundary,  i.e.,  they  give  "zero" 
reflection.  Numerical  reflections  from  the  boundaries  do  occur  for  oblique  propagating  waves; 
they  are  eliminated  by  combining  the  radiating  boundary  condition  with  increased  attenuation. 

AuemmLiga 

FOAM  treats  attenuation  forCW  propagation  by  usi.  g  a  complex  wavenumber,  k  -  (tivc)(i  ♦  t  n  P), 
where  p  =  attenuation  in  dB/wavelength,  and  n  is  a  numerical  constant  given  by  tj  =  (40  n  logjoe)*1. 
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This  is  equivalent  to  defining  an  effective  complex  wave  speed  via  i/c  — >  i/ccn  =  n/o<  1  +  .  it 
In  the  SAFE  model,  attenuation  is  implemented  for  both  comprcssional  waves  and  shear  waves  m 
a  similar  fashion:  i/cp  — >  l/cp,cfr  =  (l/cp)(l  +  i  n  Ify)  .  i/cs  — >  1/Cs.cfl  =  0/cs)(1  +  i  n  (is ) 


WITHIN  EACH  FRAME: 

-  FORWARD  PROPAGATION 

-  BACKSCATTER 

-  FULL-WAVE  SPECTRUM 
(DISCRETE  &  CONTINUOUS) 

-  SHEAR  WAVES  &  SHEAR  ATTENUATION 

-  RANGE  DEPENDENCY 

-  INCLUSIONS 


Fig.  2.  Marching-frames  approach  showing  how  SAFE  is  solved 
within  each  frame  and  the  vertical  field  from  that  frame 
-  used  as  the  starting  field  for  the  next  successive  frame. 

Marching  Frames  Approach  to  Long  Range  Propagation 

Figure  2  illustrates  how  a  FE  model  can  be  used  to  predict  long-range  ocean  acoustic  propagation. 
The  FE  full-wave  fange-dependent  Acoustic  Marching  Elements  model,  FFRAME,  (Chin-Bing 
and  Murphy,  1990)  was  the  first  of  our  FE  models  to  employ  this  technique.  A  large  "frame" 
enclosing  a  region  of  ocean  and  ocean  bottom  is  solved  by  the  FE  model.  The  vertical  column  of 
pressures  from  this  frame  is  then  used  as  a  Starting  field  for  the  next  large  FE  frame.  By 
successive  application  of  this  technique,  the  1-way  (outgoing)  propagating  field  can  be  "marched" 
to  very  long  ranges  while  still  including  the  2-way  effects  of  backscatter  losses  within  each  frame. 
Each  frame  can  be  made  a  different  size  to  accommodate  any  special  ocean  features. 
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Fig.  3.  Step-periodic  rough  bottom  with  shear  used  in  the  example. 


STEP  SIZE:  10  m  x  50  m  (WAVELENGTHS:  1/15  x  1/3) 


WATER 

1  TRAPPED  PROPAGATING  MODE  Cp=1 500  m/S 

•  SOURCE 

Density=1.0 

M  mmmmmm 

■  m  Scholte  Wave  Along  the 

|||  Cp=1800,  CS=600 

Water-Solid  Interface 

HI]  Atten,:  0.1 , 0.2  dB/w.l. 

’  ,  SOLID  BOTTOM 

Ill  Denslly=2.0 

NUMERICAL  EXAMPLE 

An  example  given  by  Schmidt  and  Jensen  (1985)  was  chosen  for  a  comparison  of  the  SAFE 
model  against  a  known  numerical  solution.  A  100  m  deep  homogeneous  water  column  (1590 
m/s)  overlies  a  solid  homogeneous  bottom  (compression  speed  =  1800  m/s,  shear  speed  =  600 
m/s,  compressional  attenuation  =  0.1  dBA,  shear  attenuation  =  0.2  dBA,  and  a  density  ratio 
between  bottom  and  water  of  2.0).  A  10  Hz.  acoustic  point  source  is  placed  at  a  depth  of  96  m 
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(-!  m  above  the  bottom)  so  as  to  excite  the  Scholte  interface  wave  that  travels  along  the  water 
solid  interface.  There  is  one  trapped  mode  in  the  water  column.  An  interference  pattern  between 
the  interface  (Scholte)  wave  and  the  one  trapped  mode  should  be  evident  in  the  waiter  near  the 
water-solid  boundary.  Roughness  was  simulated  as  shown  in  Fig.  3  using  a  step-periodic 
approach  (Evans  and  Gilbert,  1985).  Since  the  SAFE  model  can  simulate  Tange-dependent  shear, 
the  combined  effect  of  roughness  and  shear  on  water-bome  acoustic  propagation  can  be  observed. 
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Fig. 4.  SAFE  and  SAFARI  models’  TL  compared  for  a  flat  shear 
bottom.  The  deep  source  excites  a  Scholte  wave  (seen  in 
wavenumber  plot)  causing  interference  at  the  96  m  receiver. 


In  Figs.  4  and  5.  the  source  was  at  a  depth  of  96  m  (4  m  above  the  water-solid  interface)  to  excite 
the  interface  wave.  For  the  deep  receiver  (96  m)  the  interference  pattern  between  the  one 
propagating  mode  and  the  interface  wave  is  quite  noticeable  in  the  transmission  loss  (TL) 
predictions  (Fig.  4)  from  both  the  SAFE  model  and  the  SAFARI  model  (Schmidt  and  Jensen, 
1985).  A  wavenumber  spectrum  plot  confirmed  that  only  one  trapped  mode  and  one  evanescent 
(Scholte)  mode  existed.  In  Fig.  5  the  receiver  was  at  the  shallow  depth  (32  m)  and  too  far  from 
the  bottom  interface  for  the  Scholte  wave  to  interfere  with  the  one  propagating  mode.  Note  the 
absence  of  the  Scholte  in  the  wavenumber  spectrum  plot.  SAFE  and  SAFARI  agree  quite  well  for 
both  receiver  depths  (Figs.  4  and  5). 


Fig.  5.  SAFE  and  SAFARI  models'  TL  predictions  for  a  flat  shear 
bottom  -  shallow  receiver,  deep  source  (4  m  above  bottom). 


Figure  6  shows  TL  predictions  made  by  the  SAFE  model  for  a  step-periodic  rough  bottom  with 
shear  (see  Fig.  3) .  The  same  two  receiver  depths  (32  m  and  96  m)  were  used  as  in  the  flat  bottom 
with  shear  case  and  the  SAFE  TL  predictions  for  the  flat  bottom  with  shear  case  (Figs.  4  and  5) 
are  replotted  in  Fig,  6  for  comparison.  The  effects  of  both  the  Scholte  wave  and  the  step-periodic 
roughness  is  evident  in  the  TL  from  the  deep  receiver.  The  Scholte  wave  again  causes  a  "periodic 
high-frequency"  interference  with  the  single  propagating  mode  while  the  step-roughness  destroys 
this  periodic  behavior.  As  previously  noted,  the  shallow  receiver  is  too  far  from  the  bottom 
interface  to  receive  any  effect  from  the  Scholte  wave.  However,  the  step-roughness  docs  create 
backscatter  which  interferes  with  the  one  propagating  mode  resulting  in  a  "high-frequency" 
interference  effect  similar  to  that  produced  by  the  Scholte  wave  interference  effect. 
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Fig.  6.  SAFE  mode!  TL  predictions  at  two  receiver  depths,  for  two 
bottoms  (flat  and  step-periodic  rough),  each  with  shear. 


SUMMARY 

The  seismo- acoustic  FE  model.  SAFE,  has  been  developed  to  simulate  acoustic  propagation  over 
ocean  bottom  regions  that  support  compressional  and  shear  speeds  and  compressional  and  shear 
attenuations.  It  has  been  extended  to  include  marching  frames,  thus  allowing  predictions  at  longer 
ranges.  The  SAFE  model  allows  numerical  studies  of  the  effects  of  ocean  bottom  roughness  and 
compressional-shear  wave  conversions  on  the  water-borne  acoustic  pressure  field.  The  theory  and 
validations  of  the  SAFE  model  have  been  presented  and  a  numerical  example  given  that 
illustrates  the  effects  of  shear  and  interface  waves  and  roughness  on  ocean  acoustic  propagation 
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